The Drude weight of the spin-1/2 XX Z chain: 
density matrix renormahzation group versus exact diagonahzation 
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We revisit the problem of the spin Drude weight D of the integrable spin-1/2 XXZ Heisenberg 
chain using two complementrary approaches, exact diagonalization (ED) and the time-dependent 
density-matrix renormalization group (tDMRG). The aim of our paper is three-fold: First, we 
present extensive results for the temperature-dependence of D. Second, we discuss the numerical 
quality of the tDMRG data and elaborate on details of the finite-size scaling of the ED results, 
comparing calculations carried out in the canonical and grand-canonical ensembles. Third, we 
analyze the behavior of the Drude weight as the point with SU(2)-symmetric exchange is approached. 
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I. INTRODUCTION 



Due to the existence of many powerful theo- 
retical methods for one-dimensional systems such 
as bosonization^i the density matrix renormalization 
group or approaches that exploit integrabilityfi many 
properties of strongly correlated one-dimensional sys- 
tems are well studied. A paradigmatic model is the 
spin-1/2 XXZ chain which maps to interacting spinless 
fermions via the Jordan- Wigner transformation!^ While 
ground-state and thermodynamic properties (see, e.g., 
Refs. and 0) as well as many static and dynamical 
response functions (see, e.g., Refs. and H) are very 
well understood, a complete theory of transport in this 
model still poses a formidable challenge to theorists. A 
question that has attracted considerable attention for a 
long time is whether spin transport at finite tempera- 
ture is ballistic or diffusive and whether or not ballistic 
transport properties are related to integrability through 
the existence of non-trivial local conservation lawsJ^iii 
Moreover, weakly violated conservation laws can also 
cause anomalous transport properties with long relax- 
tion rates in one-dimensional systemsJ^ii^ For the most 
part, such questions have been studied in the frame- 
work of linear response theory (see Refs. [13 and [l^ and 
references therein), yet more recent work has also ad- 
dressed the non-equilibrium regime (see, e.g., Refs. [l6l - 
[2^ . Spin-transport in different experimental realisations 
of quasi one-dimensional spin systems has been probed 
indirectly by NMRS^"— and ^-SR^^ experiments. Yet, 
there is no consistent interpretation of these various ex- 
periments with regard to the question of ballistic or dif- 
fusive spin dynamics. The latest results for SrCu02, 
though, arc consistent with ballistic dynamics at elevated 
temperatures 1^ 



plementary methods, exact diagonalization and the time- 
dependent density matrix renormalization groupi^"— 
The Hamiltonian is given by: 

H^JhY. \SfSf+, + SfSf+, + A^f , (1) 



where S^, /z = x,y,z, are the components of a spin-1/2 
operator acting on site i of a chain of length L (bound- 
ary conditions will be specified later). A parametrizes 
the exchange anisotropy, Jh sets the energy scale and 
will be set to unity hereafter {h = 1). The model has 
a critical gapless phase for |A| < 1 and gapped phases 
with antiferromagnetic order and ferromagnetic order for 
A > 1 and A < — 1, respectively 1^ We will focus on the 
gapless phase |A| < 1 and study spin transport in the 
linear response regime. 

Within linear response, one decomposes the real part of 
the spin conductivity ^{uj) into a singular, zero-frequency 
part weighted with the so-called Drude weight D and a 
regular contribution at finite frequencies, CTrog(w): 



Re ct(w) = 2t:D5{uj) + CTi.cg(w) 



(2) 



Here we revisit the problem of finite-temperature spin 
transport in the spin-1/2 XXZ chain by using two com- 



Wc define transport at finite temperatures to be ballistic 
if D{T) > and our work focusscs on the dependence of 
D on temperature and exchange anisotropy A. A regu- 
lar part with a nonzero contribution to the total weight 
lo = J duj Re <j{uj) is present for A 7^ and | A| < 00. Its 
precise low-frequency behavior determines whether this 
regular channel is diffusive or not (see, e.g., Refs. [33l - l40l ). 
More recently, the momentum dependence of current re- 
sponse functions^S as well as the transverse component 
of the spin current correlations^ have also been studied. 

For the Drude weight, most studies agree that it van- 
ishes in the massive regime , ■^^^se ,42- _45 At A = 0, £> is 
finite for any temperature since the model maps onto 
non-interacting fermions. For T = 0, D can be ob- 
tained exactly via Bethe ansatz4^ Very involved gener- 
alizations of the Bethe ansatz to finite temperatures^!^ 
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consistently yield £> > for all < A < 1; however, 
the results of Refs. 113 and IH differ quantitatively re- 
garding the temperature dependence of D and the point 
of SU(2)-symmetric exchange, A = 1 (see the discus- 
sion in Refs. IH, [s^ and nSi }. Both approaches have 
also been argued to violate exact relations at infinite 
temperatures^^^ and the calculations of Ref. IH assume 
the entire spectral weight to be in the Drude peak, in- 
consistent with earlier numerical results. 



For nonzero magnetization = 
bound to D, 

< Abound < D , 



Ei^i-Sf, a lower 



(3) 



can be easily constructed from local conserved quantities 
known from the Bethe ansatz,— exploiting the existence 
of the so-called Mazur inequality. A finite bound requires 
a non-zero overlap (JQ) of the current operator J with 
at least one conserved quantity Q. The existence of such 
non-zero bounds has been used in Refs. [lol . l35l . l36l . and 
I49I to study the temperature and A-dcpendence of D at 
finite magnerizations. For S"^ = 0, however, the overlap 
of all local conserved quantities known from the Bethe 
ansatz solution with the current operator vanishes for 
symmetry reasons Only recently, a non-zero bound 
inbound > was derived by Proscn for the case of vanish- 
ing total magnetization = and < A < li^ This 
bound exploits the existence of a quasi-local conserved 
quantity that has a finite overlap with J, which can be 
viewed as the formal reason as to why D{T) is finite (see 
also Ref. [Sll for an extended discussion) . It remains un- 
clear whether or not this conserved quantity is related to 
integrability and whether the bound presented in Ref. [s^ 
is exhaustive [equality in Eq. ([3])] at all < A < 1. De- 
spite this great progress in understanding why D{T) is 
nonzero in this model, a result that previously was mostly 
based on numerical simulations , ^^'^^1"^^'^^" — there are still 
relevant open questions. First, at A = 1, there is no 
agreement as to whether or not D is finite. Second, the 
temperature dependence of D{T) for | A| < 1, albeit stud- 
ied in previous work j^^i^'^i^^i^^" — is not fully understood 
yet and is closely related to the question of which fraction 
of the total weight /g is in the Drude weight compared to 
finite-frequency contributions i^^^^i^ Finally, when eval- 
uating expressions for D with exact diagonalization, one 
can either elect to work in a grand-canonical picture (i.e., 
averaging over all subspaces with a given total 5'^) or in 
the canonical ensemble (i.e., only one sector with a fixed 
S"^ is taken into account). While one would, on general 
grounds, expect that in the thermodynamic limit both 
approaches should yield the same result, a vastly differ- 
ent behavior can emerge on finite systcms4ii^ This calls 
for a more detailed investigation. 

The first method that we employ is the time-dependent 
density matrix renormalization group for finite temper- 
atures. Various approaches for using tDMRG at finite 
temperatures have been suggested and pursued in the 
litcrature,^^"^ The main restriction of this method, 
which can otherwise provide highly accurate results, is 



the maximally accessible time scale due to the entan- 
glement encoded in the time-evolved wave-function (see 
Ref. 3 for a discussion). Here, we use the purification 
method originally proposed in Ref. |62| complemented by 
a disentangler that allows one to go to longer times than 
previously possiblci^ In a nutshell, the purification ap- 
proach relies on embedding the actual physical system 
for which one seeks to know finite-temperature properties 
into some environment. One can then show that a pure- 
state time-evolution of the full system, which is amenable 
to tDMRG, yields the thermal statistical operator for 
the system after tracing out the environment's degrees of 
freedom. For recent developments and a comparison of 
various finite-temperature approaches, see Refs. IH and 
\q^ . An advantage of this tDMRG approach is that the 
results for the accessible time-scales arc virtually in the 
thermodynamic limit. In Ref. [ssl this approach was used 
to extract the Drude weight of the spin-1/2 XX Z chain 
from the time-dependence of current-current autocorrela- 
tion functions. The results of Ref. [s^ confirm the prevail- 
ing picture that ^(r) > for |A| < 1 while for A = 1, 
the accessible time scales are too short to arrive at defi- 
nite conclusions. Here, we provide a more detailed error 
analysis of this method and we present extensive results 
for the temperature-dependence of D{T) for |A| < 1, in 
particular including A < 0. 

Our second method is exact diagonalization which 
allows to compute transport coeflacients for typically 
L <20 sitesJ^i^'^ The expression for the Drude weight 
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where Zt = X^n 6xp(— /S-Bn) is the partition function, \n) 
and En are the eigenstates and eigenvalues of H , J is the 
spin-current operator, and T = ^i{Sf Sf_^-^ + 5fS',f_|_j) 
is the kinetic energy. We present a detailed compari- 
son of a grand-canonical versus a canonical evaluation of 
Eq. A canonical evaluation of Eq. Q was system- 
atically used in a recent work^ by Herbrych et al. to 
compute the dependence of D on temperature and ex- 
change anisotropy. We argue that over a wide range of 
parameters (A and temperature T) the grand-canonical 
and canonical extrapolations yield quantitatively similar 
results. Moreover, we observe that the grand- canonical 
data often exhibits weaker finite-size dependencies than 
the canonical ones, agrees better with tDMRG data, and 
it seems to depend on system size in a more system- 
atic way throughout the parameter range of interest, i.e., 
|A| < 1. Furthermore, it has long been realized that re- 
sults for D from systems with an odd number of sites 
versus an even number of sites exhibit a different finite- 
size dependence f^i^i^ which appears to be related to 
non-trivial finite-size effects of the regular part (Treg('^) 
at small frequencies j^^j^i^ The interplay of the Drude 
weight and finite-frequency contributions was systemati- 
cally studied by Herbrych et ali^ who show that these ef- 
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fects have to be taken into account into the finite-scahng 
at incommensurate A (i.e., those values of A that can- 
not be written as A = cos(7r/i^), v integer). To avoid this 
compUcation, we focus on odd system sizes and commen- 
surate values of Aj2^ These finite-size dependencies can 
be understood by considering that the Drude weight is re- 
lated to the curvature of many-body levels upon changing 
boundary conditions via a fluxJ^iiiiH Under these three 
premises, grand-canonical evaluation of Eq. ([4]), consider- 
ation of odd L, commensurate values of A, we analyze the 
temperature dependence of D{T) throughout the gapless 
phase and we present results for D{T) extrapolated to the 
thermodynamic limit. Unfortunately, there is no theory 
for the finite-size dependence of D{T) which is the main 
drawback of the exact diagonalization approach. How- 
ever, for those parameters for which reliable results from 
ED and tDMRG are available, we report an excellent 
quantitative agreement between these two methods. To 
elucidate the behavior in the vicinity of A = 1, we follow 
D = cos(7r/z^) as a function of i' and observe that D{i') 
seems to approach finite values as -> oo (corresponding 
to A = 1). 

To summarize, the main results of our work are an 
extensive analysis of D = D{T) using both tDMRG and 
ED which are in very good agreement, a discussion of 
the numerical quality of the tDMRG data, and a detailed 
comparison of the finite-size dependencies of ED data. 

The structure of this exposition is as follows. In Sec.llll 
we summarize the expressions for the transport coeffi- 
cients that are studied in this work. In Sec. lIIIl we discuss 
details of the tDMRG approach while Sec.|lV]focusses on 
aspects specific to ED. We present our main results and 
a discussion thereof in Sec. |V] and conclude with a brief 
summary, contained in Sec. I VI I 
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FIG. 1. (Color online) DMRG data for the real-time current- 
current correlation function of an XX Z chain with A > at 
intermediate to large temperatures T. Its asymptotic value 
at long times determines the Drude weight through Eq. ((8)1. 
The calculation was carried out for a system of size L = 100 
and using a discarded weight e — 10^^ during the real time 
evolution. 



By inserting a full set of eigenstates \n) and eigenval- 
ues En in Eq. the Drude weight can be cast into 
Eq. (j4]).— 1^ This expression can be evaluated using ex- 
act diagonalization. 

The Drude weight is further related to the long-time 
asymptote of of the current-current correlation function 
{J{t)J) (see Ref.[i3): 



D = lim lim 



Re {J{t)J) 
2Lf 



(8) 



Our aim is to compute the Drude weight of the inte- 
grable XX Z Hcisenbcrg chain defined in Eq. ([Ij. The 
Kubo formula for the spin conductivity reads 



Using time-dependent DMRG, the current-current corre- 
lation function can be computed, thus yielding access to 
the Drude weight. 



1 r°° 

e-*([J(t),J])dt 



(5) 



e'^\j{t)J)dt , 



where the finite-temperature real-time current-current 
correlation function is given by 



{J{t)J)^Tv{pTe'"'je~'"'j) , pT 



(6) 



The current operator J = X); defined via a continu- 
ity equation 



i+i 



j,^tS+Sr^,/2 + R.c. . (7) 



III. DENSITY MATRIX RENORMALIZATION 
GROUP SIMULATIONS 

A. Method 

In order to evaluate Eq. ([HI) by a DMRG 
algorith m •^'^^1^^ that evolves a wave-function (as opposed 
to the time-evolution of a matrix product operator), 
one first needs to purify the thermal density matrix px 
by introducing an auxiliary Hilbert space Q such that 
Pt = TrQ|^'T)(^'T|. This is analytically possible only 
at T = oo where px factorizes. However, \^t) can be 
obtained from l^'txj) by applying an imaginary time evo- 
lution, I^t) = e--f^/(2^)|5'oo)'^'^'^ The current-current 
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FIG. 2. (Color online) The same as in Fig. [T] but for A < 
and over the full range of temperatures 0.05 < T < 5. For 
T — )■ one approaches the T = Bethe ansatz result.— The 
calculations were carried out using L — 200 as well as e = 
10"^ for T = {5, 2, 1, 0.5}, e = lO"*^ for T = {0.33, 0.2}, and 
e= 10"^ for r= {0.1,0.05}. 



correlation function of Eq. ^ can then be exactly recast 
as 

(jm = .iu/Td\ — ^ , (9) 



and this object is directly accessible in the time- 
dependent DMRG framework^^Ti^i outlined in Ref. In 
DMRG simulations, we use open boundary conditions. 

It is convenient to first express l^'oo) in terms of a 
matrix product state, 



where 



^Tcvcn ^(10) 



= (0 -1) 
= (1/V2 0)^ 



(11) 



Physical (auxiliary) sites are denoted by even (odd) 
indices. After factorizing the evolution operators 
exp(— Aif) using a fourth-order Trotter decomposition, 
they can be successively applied to Eq. ([TU|. At each 
time-step AA, two singular value decompositions are car- 
ried out to update three consecutive matrices. The ma- 
trix dimension x is dynamically increased such that at 
each time step the sum of all squared discarded singu- 
lar values is kept below a threshold value e. As men- 
tioned above, we do not only time-evolve the physical 
but also the auxiliary sites (using the physical Hamil- 
tonian but reversed time). As outlined in Ref. [ssl this 
is an exact modification to the DMRG algorithm (aux- 
iliary sites are traced over, and one can thus apply an 
arbitrary unitary transformation to them) which leads 
to a significantly slower build-up of entanglement; hence, 
longer time scales can be reached. The tDMRG approx- 
imation to {J it) J) becomes exact in the limit of e — )• 
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FIG. 3. (Color online) DMRG data for the current- 
current correlation function obtained using different discarded 
weights e during the real time evolution and different system 
sizes L. 



and AA — > 0. We ensure that both are chosen small 
enough and that the system size L is large enough to ob- 
tain results which are converged w.r.t. those parameters 
(this is discussed below). 



B. Real-time data 

tDMRG data for the current-current correlation func- 
tion are shown in Figs. [T]|4l The calculation is stopped 
once the matrix dimension has reached values of x w 
1500. This generically allows to access time scales t w 20 
which are much larger than the intrinsic scale t = 1/ Jh ~ 
1. For antiferromagnetic anisotropics < A < 1 and in- 
termediate to large temperatures T > 0.5, {J{t)J) has 
saturated for i w 20 (see Fig. [1]); thus, the Drude weight 
can be obtained via Eq. At smaller T, a large dif- 
fusive contribution with a small decay rate prevents us 
from reaching the asymptotic regimc j^^i*^^ Data for the 
isotropic chain A = 1 are shown in Figs. [T] and 01 these 
will be discussed in Sec. IV Al In the ferromagnetic case 
— 1 < A < (Fig. [2]), the current-current correlation 
function seems to saturate fast for any T. Field theory 
suggests that a diffusive contribution (of unknown to- 
tal magnitude) with an exponentially small decay rate 
existsj^i^ There are thus two possibilities to interpret 
our data at small T: either {J{t)J) saturates slowly on a 
large scale which does not manifest itself for times i < 20; 
or the magnitude of the diffusive contribution is small. 
The facts that first, we successively approach the exact 
Bethe ansatz result for T — > (see Fig. (2) and second, 
the observation of large ballistic contribution relative to 
the total weight of cr(a;) (compare Sec. IV B|) support the 
latter, and in the following we will rely on this interpre- 
tation and discuss the available numerical data. 

As mentioned above, we have to ensure that the dis- 
carded weight and the Trotter step size are chosen small 
enough and that the system size L is chosen large enough 
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FIG. 4. (Color online) The same error analysis as in Fig.[3]but 
for A = ±1. Note that for T = oo curves at ±A conincide.— 



to obtain results which are converged w.r.t. those pa- 
rameters. Typical values are L = 200, AA = 0.2, and 
e — 10"^ as well as e = 10"'* . . . 10"^ during the imagi- 
nary and real time evolutions, respectively. The depen- 
dencies on L and e are illustrated in Figs. [3] and HI The 
finite system size manifests itself on time scales t > L/v 
with V being a proportionality factor which has the di- 
mension of a velocity. This is the so-called light-cone 
effect in nonrelativistic systems: correlation functions 
{A{t)B) of local operators A and B fall off exponentially 
for X > vt with x denoting the spatial distance between 
the regions on which A and B act.— Note that in our case 
for times t <20 and temperatures T > 0.5, the data for 
all L > 75 coincides (see Fig. Moreover, if L is cho- 
sen too small or if e is chosen too large, one generically 
underestimates the Drude weight. 



IV. EXACT DIAGONALIZATION 

Our ED analysis of the Drude weight is based on an 
evaluation of Eq. (j?]). We use periodic boundary condi- 
tions and we compute transport coefficients for systems 
as large as L < 19. We will make a distinction between a 
grand- canonical versus a canonical evaluation of Eq. ((4|), 
which includes the evaluation of the kinetic energy from 



(12) 



To be specific about grand- versus canonical interpreta- 
tion, the sums over eigenstates in Eqs. ^ and have 
to be read as: 
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FIG. 5. (Color online) ED data for the leading coefficient 
C{L) of D(T) « C/T in the high-temperature limit (see 
Eq. (|13p ). (a) Grand-canonical and (b) canonical data at 
A = 0.5, comparing data from odd L to even L. The lines are 
an extrapolation in system size using C{L) = C + A/ L + B / L'^ 
using data-points with 5 < L < 19. (c) Comparison be- 
tween grand-canonical and canonical data for C for odd L 
and A = 0, 0.5, 1. The ED results are extrapolated to L = oo 
using Eq. ([T4)l . 



in the former case, whereas in the latter case, we fix 
to a desired value and sum over all eigenstates En{S^) in 
that subspacc only. For even system sizes and in the 5*^ ~ 
subspacc, the Drude weight vanishes on finite systems. 
Therefore, the smallest that can be considered in a 
canonical calculation is S*" = 1/2 on odd system sizes 

Let us first compare data obtained from chains with 
an even versus an odd number of sites in the example of 
A = 0.5 and in the high-temperature limit. We compute 
the leading coefficient of D{T) = C/T + OiT''^) in a 
high-temperature expansion in 1 /T from^ 



C - hm [TD{T)] 

T-i-oo 



(13) 



The results are shown in Fig.[5l^a) for the grand-canonical 
case and in Fig. [Hb) for the canonical one. In the latter 
case, odd system sizes correspond to working with = 
1/2 whereas for even system sizes, we work in the = 1 
subspace (note that S'^/X — )• as L — >• oo)»^ The figure 
further contains fits to the data using 



C{L) + A/L + B/L^ 



(14) 



The results for C are all consistent with each other, with 
a somewhat larger deviation for the canonical data with 
= 1. However, the coefficient B is larger for chains 
with an even number of sites, which is particularly evi- 
dent in the grand-canonical case [see Fig.lSJa)], indicating 
that the data for odd system sizes approaches the large-L 
limit faster at T — > cx). This behavior extends down to 
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FIG. 6. (Color online) ED data for the Drude weight D 
for A — 0.5. The filled (open) symbols are the results for 
the grand-canonical (canonical) calculation of D. The solid 
(dashed) lines are an extrapolation of D in system-size using a 
second-order polynomial in 1/L, Eq. (|15p. for grand-canonical 
(canonical) data. 



temperatures of T ^ 0.5. In the following, we therefore 
restrict our analysis to odd L = 5, 7, . . . , 19. 

Next we compare grand-canonical data to canonical 
data for C for A = 0, 0.5, 1 which is shown in Fig. [S^c), 
extrapolated using Eq. (jH)) . We observe that in the 
high-tempcrature limit, B w 04^ Generally, one would 
expect both the grand-canonical and the canonical data 
sets to result in the exact same number as i — > oo. 
This is consistent with our data for A = 0. Clearly, 
though, the grand-canonical data exhibts no finite-size 
dependence at all for A = in contrast to the canonical 
data. Therefore, we conclude that data obtained from 
a canonical evaluation of Eq. approaches the large L 
limit slower than grand-canonical data. At A = 0.5, the 
grand-canonical data and the canonical data extrapolate 
to the same value. At A = 1, the grand-canonical data 
still follows the simple function given in Eq. ([H)) quite 
well and extrapolates to a finite value of C'{L ~ oo) (this 
reproduces the results of Refs. [l^ and [13), whereas the 
canonical data do not follow any low-order polynomial 
in 1/L. The authors of Ref. [H^ argue that the canoni- 
cal data for A = 1 can only be interpreted in terms of 
C{L = oo) = 0. Our point of view is different: since there 
is no theory for the finite-size dependence of D, one can 
only use the canonical data at A = 1 extracted from the 
largest L as an upper bound to the true D, assuming that 
the finite-size dependence is monotoneous. To summa- 
rize, while the grand-canonical data can be well fitted to 
a simple polynomial in 1/L of low order for all |A| < 1, 
it is at present unclear how to extrapolate the canonical 
data close to A = 1 to L -> oo. 

The main goal of our work is to obtain D = D{T) 
from exact diagonalization by means of an extrapolation 
in in order to compare with tDMRG data. In Fig. [HI 
we demonstrate, for the example of A = 0.5, that the 
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FIG. 7. (Color online) Comparison of tDMRG and ED data 
(extrapolated in L) for the Drude weight for (a) A = 0.5, (b) 
A = cos(7r/4), and (c) A = cos(7r/5). Circles are the exact 
T — Q Bethe ansatz results taken from R ef . l4a . the stars are 
the T > Bethe ansatz results from Ref. Be- 



fitting function 

D{T,L) ^ D{T) + A{T)/L + B{T)/L^ (15) 

describes the available data relatively well at all tempera- 
tures, which has previously not been appreciated in exact 
diagonalization studies (Herbrych et al. also present re- 
sults for D(T) from such extrapolations as a function of 
T in Ref. l54l ). However, in the temperature regime of 
0.1 < T < 0.5, the monotony behavior of D{T,L) as a 
function of L changes dramatically, which is particularly 
evident by comparing data from odd to even L (see the 
discussion and figures in Ref.Ei, rendering this regime 
the most difficult one and the least reliable for finite-size 
extrapolations. In Ref. [s^. it was observed that D = C/T 
approximates the full D{T) well down to temperatures of 
T ^ 0.5, suggesting that it should be possible to obtain 
an estimate for D{T) from finite-size extrapolations in 
the regime 0.5 < T < oo. As T 0, we observe A ^ 
(see the discussion in Ref. IH and references therein). In 
the following, we will use Eq. (fTSjl to obtain D{T) for 
several A = cos(7r/i/), with z/ = 3, 4, . . . , 25 and at A = 1 
(corresponding to v ^ oo). Error bars in the figures 
are obtained from excluding either the largest or small- 
est system size from the fits. We have also tried other 
possible extrapolation functions, yet we find that a sec- 
ond order polynomial in 1/L consistently approximates 
the data the best. 



V. RESULTS AND COMPARISON OF TDMRG 
WITH EXACT DIAGONALIZATION 

In this section, we present our main results, namely 
ED and tDMRG data for the temperature dependence of 
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FIG. 8. (Color online) Comparison between tDMRG and ED 
data for A — —0.5. 
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D{T), a comparison of canonical versus grand-canonical 
ED data, and wc revisit the case of A = 1, for which no 
exact results are available on whether D{T) is finite or 
not. Note that results for D{T) from an extrapolation of 
canonical ED data were previously presented in Ref. [H 
(we quantitatively agree with their data for A = 0.5), 
while grand-canonical ED data for finite systems was ex- 
tensively discussed in Ref. IH (see also further references 
mentioned in that work). 



A. Gapless regime: |A[ < 1 

Figure [7] contains representative examples of our re- 
sults for D{T) for three positive values of A = cos(7r/i') 
{v = 3, 4, 5). We emphasize the main observations. First, 
the canonical and grand-canonical ED data agree well 
with each other both at low and high temperatures. At 
low temperatures, the subspacc with the smallest value 
of S'^ dominates the grand-canonical sum in Eq. ^ and 
since this smallest S'^-subspace enters in the canonical 
calculation of D, the agreement at low temperatures is to 
be expected. Second, the deviations between the canon- 
ical and grand-canonical data are the largest at interme- 
diate temperatures T ^ 0.5, and increase as v increases. 
Third, the agreement between the tDMRG data and the 
ED data is excellent, and the tDMRG appears to be 
closer to the grand-canonical ED curves. In all cases, 
the extrapolated ED data are consistent with the exactly 
known zero-temperature result^ (circles in Fig. [7]). 

Our grand-canonical data further agree with the avail- 
able finitc-tempcraturc Bcthc ansatz results from Ref. \^ 
at T > 0.4 (while the canonical data exhibit larger devia- 
tions). Yet, upon lowering temperature, the Bethe ansatz 
results are systematically below the ED data and follow 
a quantitatively different temperature dependence with a 
non-zero slope at T = (see also Ref. |5a). Our ED data 
suggest a vanishing slope of D(T) as T — >■ 0, in agreement 
with the results of a different finite-temperature Bethe 
ansatz calculation4^ This discrepancy can at present not 
be resolved yet may be amenable to alternative tDMRG 



FIG. 9. (Color online) (a) Drude weight versus v, where A = 
cos (^). (a), (b) Comparison of tDMRG (filled symbols) and 
ED (open symbols) using grand-canonical data in a log-log 
plot. Thin solid lines connecting the ED data are guides to the 
eyes, (b) C vs u {C ^ liniT-^oo [T -D(T)]). The thick dashed 
line in (a) shows an extrapolation using D{u) = A ■ (i^)^^ 
using v = 3, . . . ,7. 



methods for finite temperature.— >^i^ Note, though, 
that if the Drude weight at A = 1 is finite, then its 
zero-temperature slope dD/dT has been predicted to be 

As an example for the behavior of D{T) in the fer- 
romagnetic, critical regime — 1 < A < 0, we present 
results for A = —0.5 in Fig. El In this case, we were 
able to extract accurate tDMRG data down to very low 
temperatures (compare the discussion of Sec. lIIII and the 
real-time data presented in Fig. The agreement be- 
tween ED data on the one hand - both canonical and 
grand-canonical data - and tDMRG on the other hand is 
excellent, at all temperatures considered. 

We next study how D{T) behaves as A = 1 is ap- 
proached, sending v to large values. We compute D{T, v) 
from grand-canonical ED data for various values of i/ 
from an extrapolation of D{T, v, L) in system size that is 
carried out for each pair ofT^v, along the lines discussed 
in Sec. |IV] (see Fig. [5]). Thus, for each temperature T, 
we obtain a series of data points D(T, v) as a function 
of I/, where A — >■ 1 as i/ — >■ oo. Examples are shown in 
Fig. inja). For v < 10, we include tDMRG data that is 
in reasonable qualitative agreement with the ED results 
(for 1/ > 10, we cannot reliably extract results for D{T) 
from the real-time tDMRG simulations). The figure fur- 
ther shows that the data do not follow a power-law in h' 
(dashed line in the figure), but rather, D{T,v) seems to 
level off at finite and infinite temperature, suggesting a 
non-zero value of D(r) at A = 1. The real-time DMRG 
data for A = ±1 (Figs.[T]and|3]) do not contradict a possi- 
bly finite Drude weight; with the availiable time scales it 
is not possible to discern between the case of D > or a 
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FIG. 10. (Color online) Fraction of the Drude weight 
2itD/Io to the total spectral weight of the spin conductiv- 
ity, 7o, obtained from ED by grand-canonical calculation of 
Eq. ((4|. lo = J Re a{uj)duj, through the f-sum rule given by^ 
lo = j;{—T). The extrapolation to 1/ — >■ oo is done as in 
Fig. El 



current correlation function that decays on a fairly large 
time scale which cannot be accessed by our numerics. 



B. Spectral weight of the Drude peak 

Finally, we study how much of the total weight Iq of 
the spin conductivity is in the Drude weight, compared to 
finite-frequency contributions. A similar analysis of ED 
data was presented in Ref. [13 using canonical ED data 
at T = oo exploring the dependence on A and , which 
we here complement by using grand-canonical data and 
consideration of finite temperatures < T < 2. 

The real part of the spin conductivity obeys a sum 
rule: 



Iq 



duj Rea{uj) = — ( 
L 



-f). 



(16) 



We plot 27rZ)//o in Fig. [TO]as a function of temperature 
for A = —0.5,0.5 and 1. These data are obtained from 
grand-canonical ED results plus an extrapolation using 
Eq. ((T5|) . The units are chosen such that the vertical axis 
can be read in %, implying that 2ttD/Iq = 1 corresponds 
to the absence of any finite-frequency contributions. 

At A — ±0.5 and in the low temperature limit, practi- 
cally all the weight sits in the Drude weight which then 
monotoneuosly drops to values larger than 50% of the 
total weight at high temperatures (at T = oo, the values 
for A = 0.5 and —0.5 should agree with each other).— 
By all means, the contribution of the Drude weight is 
therefore substantial at all temperatures. At A = 1, the 
low-temperature regime is difficult to obtain from ED 
data due to several reasons. It is well-known that the 
low-temperature behavior of thermodynamic quantities 
at A = 1 is strongly infiuenced by the appearence of 



a marginally relevant perturbation to the effective low- 
energy field theory, which otherwise, in the gapless phase, 
is simply a Luttinger liquid. This cannot be captured by 
ED on short chains (compare the work by Lafiorencie et 
ai, Ref. ItiI ). and we thus expect the temperature de- 
pendence of the true Drude weight to deviate from our 
result at low temperatures T < 0.1. The exact T = 
result for 2hD/Iq (circle) sits well below the ED data4^ 
Based on the available ED data, we conjecture that one 
possible behavior of D{T) and 2Tr D/Iq in the thermody- 
namic limit at A = 1 is a maximum at some intermediate 
T ~ 0.1 if D(T) is non-zero at all. This would imply a 
positive slope of D(T) at zero temperature, in agreement 
with Refs. [43 and l56|. The existence of the additional 
maximum in D(T) ai T 0.8 cannot be resolved with 
the available system sizes and should thus be considered 
as spurious. An intriguing observation from Fig. [10] for 
A = 1 is that the Drude weight amounts to 80% of the 
total weight at temperatures T ^ 0.1, which drops to less 
than 10% at high temperatures. 

Our results indicate that the contribution of the Drude 
peak is substantial for all A > and A sufficiently 
smaller than 1, and may even be large at low temper- 
atures (see the Fig. [10] with A = 0.5). We further sug- 
gest that our data for the temperature dependence of 
2ttD{T)/Iq at a = ±0.5 could be used to obtain esti- 
mates on the diffusion rate 7 introduced in Refs. Iss] and 
[36l along the lines of Ref. or for a comparison with 
QMC results 1^ Finally, we observe that for the commen- 
surate values of A = cos{tt/v) studied in our work, the 
lower bound to D{T) from Prosen's work^ seems to be 
almost exhaustive. 



VI. SUMMARY 

In this paper we studied the problem of the spin Drude 
weight of the spin-1/2 XXZ Heisenberg chain at finite 
temperature using the time-dependent density matrix 
renormalization group as well as exact diagonalization. 
This complements earlier works in various ways. First, 
we elaborated on numerical details of both methods. 
tDMRG yields the real-time current-current correlation 
function ( J(t) J), and D can be extracted from its asymp- 
tote provided that the asymptotic time regime where 
{J{t)J) saturates can be reached. We presented data 
for various discarded weights which controls the accu- 
racy of tDMRG calculations. While the thermodynamic 
limit L — > 00 can be accessed easily within tDMRG, 
exact diagonalization is bound to system sizes L < 20 
and thus requires an extrapolation to L — > 00. We 
showed that canonical and grand-canonical ED data are 
in good agreement at small A < 1. Grand-canonical 
data appear to have more systematic finite-size depen- 
dencies all across the gapless phase. As a second key 
aspect of our work, we presented extensive data for the 
temperature-dependence of the Drude weight. For pa- 
rameters where the asymptote of {J{t)J) can be accessed 
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by tDMRG (antiferromagnetic chains at intermediate-to- 
high temperatures as well as for ferromagnetic chains at 
any T) , results are in perfect agreement with ED extrap- 
olated to L — >■ c». Finally, we investigated how D{T) 
behaves in the vicinty of A = 1 by sending v to large 
values (A = cos(7r/z^)). These data are consistent with 
p(T) > at A = 1 (in disagreement with Refs. 113 and 
|54[ ). and a final clarification of this open question is left 
for future work, exploiting, for instance, alternative and 



novel finite-temperature tDMRG mcthods i^^'^'^i 
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